Tree-Level Electron-Photon Interactions in Graphene 
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Graphene's low-energy electronic excitations obey a 2+1 dimensional Dirac Hamiltonian. After 
extending this Hamiltonian to include interactions with a quantized electromagnetic field, we calcu- 
late the amplitude associated with the simplest, tree-level Feynman diagram: the vertex connecting 
a photon with two electrons. This amplitude leads to analytic expressions for the 3D angular depen- 
dence of photon emission, the photon-mediated electron-hole recombination rate, and corrections to 
graphene's opacity 7ra and dynamic conductivity 7re^/2/i for situations away from thermal equilib- 
rium, as would occur in a graphene laser. We find that Ohmic dissipation in perfect graphene can 
be attributed to spontaneous emission. 
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Electron-photon interactions determine the opto- 
electronic properties of a material. The electrons in 
graphene, a single atomic layer of graphite, exhibit su- 
perlative electronic properties associated with their ex- 
otic Hamiltonian [l!, 01- In particular, a tight binding 
model Q of graphene produces a Hamiltonian that, for 
low energy excitations, is formally identical to a 2-1-1 di- 
mensional Dirac equation for massless fermions with 
the Fermi velocity and the sublattice state vector filling 
the roles of the speed of light and spin respectively. As 
part of an effort to understand how electron-hole recom- 
bination might limit the function of a graphene-based 
transistor, we use this Dirac Hamiltonian to calculate 
the amplitude for the electron-photon interaction dia- 
grammed in Fig. [TJ Rotating this diagram with respect 
to the time axis allows the consideration of both pho- 
ton emission (i.e. recombination) and absorption rates, 
which we relate to graphene's opacity and dynamic con- 
ductivity. 

These measurable 043 properties have been previ- 
ously treated using semiclassical methods (where the 
electromagnetic field is not quantized) within the Knbo 
and Landauer formalisms |9l4l3l[ and perturbation the- 
ory Our fully quantum mechanical calculation 
reproduces results found previously, such as ira for the 
optical opacity [1,H, 15| and ire^ /2h [lol- 15| for the zero- 
temperature conductivity. We extend these previous re- 
sults to non-eqnilibrium situations {e.g. population inver- 
sion) and specify the full angular dependence of photon 
emission/absorption. Furthermore, we identify sponta- 
neous emission as the mechanism of dissipation, present 
even in idealized graphene, that is usually left unspecified 

The carbon atoms in graphene form a two-dimensional 
honeycomb network with two ineqnivalent atomic sites 
per unit cell. In the simplest tight-binding description of 
graphene, an electronic energy £ is associated with each 
atomic site in the sheet, and an energy t parametrizes 




FIG. 1: Schematic drawing of a representative emission pro- 
cess (left), and the corresponding Feynman diagram (right). 
The photon lives in 3D space, while the electrons are confined 
to the graphene sheet. The initial electron is described by its 
momentum Pi = ftki and its pseudospin =>, which for a con- 
duction electron near K"*" is directed along k^. Interacting 
with the photon (wavevector k^ and polarization e) destroys 
the conduction electron, creating a valence band electron with 
momentum /ik/ and pseudospin <=. 



the probability of an electron hopping from one site to 
its neighbor on the other sublattice. An operator 
creates a 2Pz electron on the 'A' site in cell j, with a 
corresponding destruction operator . . With similar 
operators for the 'B' sites, the total Hamiltonian H is 

H = £Y. (^R, ^R, + ^Rj ) - * E (^Ri ^Rj + ^•'^•) , 

i <i,j> 

(1) 

where j runs over the N sites in the sheet, and i runs over 
the nearest neighbors of the site j. Spin indices on the op- 
erators and the sums are understood. Fourier transform- 
ing the creation and annihilation operators {e.g. A^i- — 
J2j ^Q, exp(iR,-Qj)/y]V, where the Qj = ^bi + ^^^ba 
are the N — N1N2 wavevectors in the first Brillouin zone) 
allows the Hamiltonian ([1]) to be written. 
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FIG. 2: The hexagonal first Brillouin zone (left) and the dis- 
persion relation near the points (right). On the left, the 
K"*" points are indicated by thin arrows, and the reciprocal 
lattice primitive vectors bi by thick arrows. Shading indi- 
cates how translating some slices of the hexagon by reciprocal 
lattice vectors reconstructs an equivalent Brillouin zone, here 
shown in a bowtie configuration, that centers the inequivalent 
points in two triangular regions. Near the points the 
dispersion relation is linear in |k|, which gives the Dirac cones 
shown on the right. Absorption or emission of a photon trans- 
fers an electron from one cone to the other. 



There are two spin states per Qj, and two mobile 2Pz 
electrons per cell, so the first Brillouin zone (Fig. [5]) is 
exactly filled in electrically neutral graphene at zero tem- 
perature. The energy origin is set at the energy of the 
highest occuppied states, which are those at the Brillouin 
zone corners K'' = k^^2_!±i + ^bi -h nbs [l3|. The label 
K = ±1 indexes the two inequivalent corners. For Q near 
a K'' point the single-particle Hamiltonian T-L is 



T-L = VpiKCTxPx + (^yPy), 



where the momentum p = ?ik = /i(Q — K). With 
defined by the in-plane components p of p = p -\- PzZ 
according to p = p(cos(/)x -|- sin</)y), the corresponding 
eigenvalue equation is 
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(4) 

where the band index /3 = ±1 labels whether the en- 
ergy is positive or negative {i.e. conduction or valence). 
Thus the Hamiltonian ([3]) produces a linear dispersion 
relation £ = ±vfP. The product /3k/2 gives the helicity 
eigenvalue for the state |x), where the helicity operator 
is defined as h — 'H/{2hvfp). 

Having identified the eigenspinors |x(p, k)) of the 
single particle Hamiltonian we can re- write the total 
Hamiltonian H in terms of operators that create (Cq) 



and destroy (Cq) energy eigenstates, 

H = Y,^VFp{{Xc\Clci + {Xv\Clci) 
Q 
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(5) 



{\Xc)Cc.Q + \Xv)a.Q) 



where the sum is over Q near K'' and c (v) refers to the 
conduction (valence) band. 

We introduce the electromagnetic field with a Peierls 
(minimal coupling) substitutio n p — p — qA/c, treating 
the new vector potential term [17| as a quantized pertur- 
bation H' in the full Hamiltonian H = Hq + H' , 
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Here j indexes the photon's polarization states, V is the 
normahzation volume, is the relative permittivity, and 
uj — c|k|. As is evident from the appearance of the speed 
of light c (and not the Fermi velocity vp) in this sub- 
stitution, the electron-photon coupling implied follows 
from the local gauge invariance of the standard model 
Lagrangian, and is not related to the properties of the 
Hamiltonian ([3]) under gauge transformations. 

The electron-photon interaction rate can be calculated 
using the standard arguments of Fermi's Golden Rule, 
suitably modified to account for the system's mixed di- 
mensionality. The rate Ti^j to go from an eigenstate 
\ipi) of the unperturbed electronic Hamiltonian Hq to a 
given final state \(pf) is 



r,^y--|(^/(t)|^(t))|2, where 



(7) 
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{^fit')\n'{t')\ip,{t'))dt'^M. (8) 



In the position representation, the time-dependent solu- 
tions to the unperturbed electronic H have the form 

(rX<)) - -^e^^'^-''-^'^F{z)\x{hk,f3,n)), (9) 



where lo = VF\k\ and A is the graphene area. Initially we 
consider processes that create a valence electron \ipf) oc 
\x{fiky, —1, k)) and a photon |k^,e), while destroying a 
conduction electron \ipi) cx \x{hkc, 1, k)). Then 
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Derived from the z-extent of the carbon 2Pz atomic or- 
bitals, the normalized function F{z) is only appreciable 
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within a few angstroms of the graphene plane. Since 
we are considering photons with optical or longer wave- 
lengths A — 27r/fc.^, the integral over dz' gives unity to 
excellent approximation. In atomic physics this step ap- 
plies to all three spatial dimensions (e*'' "" ~ 1) and is 
known as the dipole approximation. 

We square M, and consider the interval Ai to be short 
compared to the lifetime l/F and long compared to the 
time scale 1/w set by the energy of the transition, i.e. 
1/r 3> At ^ In this limit, with large area A, 
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where we have used the standard relations |n) — 
^/l±n\l±n), C\n) = ^/n\±{n~l)), and {n'\n) = Sn'n, 
with the upper (lower) sign chosen for bosons (fermions) . 
Thus we see that the recombination rate is proportional 
to the number of conduction electrons ric and the num- 
ber of holes 1 — riv The first and second parts of 1 -I- n-y 
correspond to spontaneous and stimulated emission re- 
spectively [18 1 . 

To evaluate the angular matrix element in 
(fTTj) . we define an orthonormal triple k-^ — 



(sin 9^ 



C0S( 



, sm f T, sm ( 



, cosf 



y), £i = z X k-., 



and €2 — k^y. £i/|k^ x e^l that describes the photon and 
its possible polarizations. Summing over the possible 
polarizations j of the created photon gives 

Y.\{Xv\c}-i*\Xc)? = 1- sin2 sin2 (0,/2 + 0, /2 - ) . 

(12) 

As Pz does not appear in %, the component of the photon 
polarization along z does not contribute to this matrix 
element. The integrals over the energy and momentum 
(5-functions in (jlip can now be performed, with the result 
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The last line in (jl3p corresponds to the angular matrix 
element p^ . 

Since vp/c is a small number ^ 1/300 Q, several ap- 
proximations are in order. To better than 1% accuracy 
kc ~ ky and k-y/kc(v) — 2vf/c. The energy of the ini- 
tial conduction electron is half that of the photon, and 
of the same magnitude but opposite sign of the final va- 
lence electron. The photon's momentum is negligible in 




FIG. 3: Polar plot of photon emission distributions in the xz 
plane. The dashed (dotted) curve corresponds to the emis- 
sion from an initial electron moving along the x-axis (y-axis), 
while the solid black curve represents the average over all ki 
directions. The Lambertian function cos 6 is shown in grey for 
comparison. The inset shows the 3D pattern for one choice of 
k. An electron moving along y emits a x-polarized photon, 
and thus cannot emit in the x direction. 



comparison to the electrons'; as a result (jic ~ and 
K+ o transitions are impossible in this low energy 
limit. To lowest order in vp/c the angular dependence of 



(Uni) is (l-sin^e*^ sin^ 



- 0c) - 2^ sin 0^ cos(0^ -(kc))- 



Thus for small vp /ca. conduction electron is slightly more 
likely to emit a photon opposite kc than along kc- Fig- 
ure [3] shows various plots of the angular distribution in 
the small vf/c limit, which we will adopt henceforth. 

When averaged over the possible momentum directions 
of the conduction electron, the emission or absorption of 
a photon depends on the polar angle 9~f from the normal 
to the graphene sheet like 1— ^ sin^ 9^. Because this func- 
tion falls off more slowly than the Lambertian function 
cos 9, a graphene sheet will appear progressively brighter 
{i.e. blacker) at angles away from normal incidence. At 
this level of analysis the angular matrix element (|12l) . 
and thus the rate, is zero for the metallic nanotube case 



20| . The interaction Hamiltonian contains only pho- 
tons polarized along the nanotube axis, and such photons 
do not couple the initial and final electronic states. 

The form of the matrix element ((T^ indicates that an- 
gular momentum conservation is enforced in an unusual 
way in this problem. In a more conventional condensed 
matter system a typical optical transition involves bands 
with different orbital angular momentum quantum num- 
bers, and allows the possibility of a spin flip. For in- 
stance, in gallium arsenide interband transitions occur 



between orbitals with S and P symmetries [2l|. Here 
the transition is 2Pz — 2Pz and there is no spin flip. 
Thus the usual sources of angular momentum for the 
photon do not contribute in graphene. The structure 
of the matrix element (|12|) . which follows directly from 
the Hamiltonian ([3]) and the assumption of minimal cou- 
pling, implies that the pseudospin flip creates the angu- 
lar momentum h of the photon. We further explore this 
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connection between pseudospin and angular momentum 
elsewhere po| . 

For states connected by the (5- functions in Eq. ([TT|) , the 
proportionality T^^y oc {1 — Uy) [1 + rij) is general, 
and applies whether the n's reflect equilibrium distribu- 
tions or not. Non-thermal distributions are commonly 
handled by introducing a quasi-Fermi level that differs 
for electrons and holes [22j . To simply illustrate the time 
scales relevant for photon-mediated electron-hole recom- 
bination, we consider a perfect population inversion, i.e. 
He = 1 and n„ = 0. Integrating (fT3|) over all directions of 
gives the rate for a conduction electron with energy 
£c = ftWc = vpfikc to decay spontaneously (n^ = 0) via 
photon emission. 



stimulated emission rate cx ric (1 — n„) n^, which gives 



8a 
3er 



(14) 



where a — /he ~ 1/137 is the fine structure constant. 
This rate corresponds to a lifetime r = l/F of about 3 ns 
for a 1 eV conduction electron. 

For thermal populations the averaged transition ma- 
trix element (|Mp) = TY{pcPvP-i\M\'^}, where the den- 



sity operators p are given by p 



-n/kT 



/Tr{e 



-n/kT 
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and the trace is taken over the possible occupations: 
n S [0, 1] for the electrons, and n £ [0, oo] for the photons 
p3| . Evaluating the trace gives Bose and Fermi distribu- 
tion functions, 
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(15) 



where we have allowed for a chemical potential p. The 
second line of (fT5|) shows that recombination stimulated 
by the blackbody background becomes important for 
2£c = ^ ^ kT. At room temperature with p — a 
conduction state with energy £c — kT ~ 0.025 eV will 
be populated and decay with a characteristic lifetime of 
about 400 ns. For many practical purposes this rate is 
negligible, since, for instance, the second order (Auger) 
process gives picosecond lifetimes (2^ . 

In contrast, the time-reverse of this recombination pro- 
cess, photon absorption, is observable practically by the 
unaided eye ^, ^8]. To analyze absorption we proceed 
as in the derivation of Eqs. (fT3l - [T4| . this time consider- 
ing illumination normally incident on the graphene plane 
{6j = 0) at a rate = njcA/erV ■ Then the promotion 
rate r^_>c from the valence to the conduction band is 



Ty^c = Tra (1 - He) nyTi, 



(16) 



where we have included a factor of 4 for the valley and 
(normal) spin degeneracies. Discounting spontaneous 
emission into the illuminating beam, we take the net ab- 
sorption rate Tabs to be the promotion rate minus the 



(17) 



With initial ric = and n„ = 1 Eq. p7|) reproduces the 
ira result for the optical absorption of a graphene sheet 
0, 0] , and identifies spontaneous emission as the source 
of dissipation. For ric > Uy the absorption is negative, 
implying gain and the possibility of a graphene laser 18, 
[25| . As before, thermally averaging pT| replaces the n's 
with Fermi functions, with the result that the absorption 
goes to zero for fku ^ kT or huj <^ 2|/i|. 

We can relate the energy absorption rate implied by 
P?)l to the conductivity a by invoking Ohm's Law, 
which implies that the power dissipated per unit area 
is K • E = (tE^. Here K is the current density and 
E = —^dA/dt is the electric field. Since the energy den- 
sity of the electromagnetic field is erE^/47r = n^hio/V, 
we have 



cue Tre 
(7 = — (ny - n-c) ^ — {uy - Uc), (18) 



which is 7re^/2/i at T = 0. This expression can be written 



Ah 



, , huj + 2u\ f hw — 2ii 

tanh + tanh ' 



4fcT 



V AkT 



(19) 



after thermal averaging, which is identical to the result 
found previously 0,M ■ Our calculation, like the pre- 
vious ones, is not rigorous at w = 0, as the dc limit ex- 
plicitly violates the assumption required to generate the 
energy (5- function in Eq. (|lip . 

In conclusion, we have performed the first calculation 
of graphene's optical properties with a quantized electro- 
magnetic field. The calculation is fully quantum mechan- 
ical and free of thermodynamic assumptions until the fi- 
nal step, which allows the treatment of systems far from 
thermal equilibrium. Furthermore, the inherently three- 
dimensional formalism gives amplitudes as a function of 
photon polarization and propagation direction relative 
to the graphene plane. The dependence on electron and 
photon state occupation numbers follows directly from 
their fermionic and bosonic commutation relations. The 
former result in Pauli blocking, while the latter give terms 
that can be identified with spontaneous and stimulated 
emission. Spontaneous emission proves to be a source of 
dissipation present even in idealized graphene, with an 
implied violation of time-reversal symmetry whose intro- 
duction can be traced back to the use of Fermi's Golden 
Rule. Stimulated emission from graphene could prove 
technologically useful, since an electronic population in- 
version would allow graphene to serve as the gain medium 
in a laser tunable over a broad band of frequencies. 
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